DESY 98-132 
HUB-EP-98/60 
TPR-98-28 
September 1998 

1 



00 
On 

00 
(N 



> 
00 

t> 

OS 

o 

oo 

D 



X 



Renormalization of four-fermion operators for higher twist calculations* 

S. Capitani a , M. G6ckeler b , R. Horsley c , H. Perlt d , P. Rakow b , G. Schierholz a <° and A. Schiller d 
a Dcutsches Elektronen-Synchrotron DESY, Notkestrasse 85, D-22607 Hamburg, Germany 
b Universitat Regensburg, Institut fur Theoretische Physik, D-93040 Regensburg, Germany 
c Humboldt-Universitat, Institut fur Physik, Invalidenstrasse 110, D-10115 Berlin, Germany 
d Universitat Leipzig, Institut fur Theoretische Physik, Augustusplatz 10/11, D-04109 Leipzig, Germany 
c Deutsches Elektronen-Synchrotron DESY and NIC, Platanenallee 6, D-15738 Zeuthen, Germany 

The evaluation of the higher twist contributions to Deep Inelastic Scattering amplitudes involves a non trivial 
choice of operator bases for the higher orders of the OPE expansion of the two hadronic currents. In this talk we 
discuss the perturbative renormalization of the four-fermion operators that appear in the above bases. 



1. INTRODUCTION 

We present here the first steps of a more general 
and large scale research program, which aims at 
the evaluation of the contributions of higher twist 
operators in Deep Inelastic Scattering (DIS) pro- 
cesses |l|-^]. We consider in particular the next- 
to-leading terms in the OPE, which are the next 
order corrections to the presently known DIS am- 
plitudes and structure functions. The matrix el- 
ements of higher twist operators are outside the 
reach of perturbative QCD, and lattice QCD pro- 
vides at present the only way to compute them in 
a reliable way and from first principles. Renor- 
malization factors are then necessary to relate the 
numbers extracted directly from the lattice to the 
physical matrix elements. 

2. HIGHER TWIST 

In the light-cone expansion of the T-product 
of two hadronic currents in DIS, the leading be- 
havior of the Wilson coefficients is given (up to 
logs) by C n ^(x 2 ) ~ ( x 2)(d o( „, l) -™-2d J )/2 ; and ig 

governed by the twist r = d,Q{ n ,i) — n (dimension 
minus spin) of the corresponding operators. The 
various terms in the expansion can thus be classi- 
fied according to their twist: operators of twist-2 



are the leading ones, and the higher twist contri- 
butions are suppressed^ as (gf) <T T """ )/2 . 

Our main interest is A T for the next-to-lowest 
twist operators. Twist-4 operators of spin n are 
related to the 1 /Q 2 power corrections to the n-th 
moment of the structure functions: 



dxx- 2 F 2 (x,Q 2 ) = C( 2 )(g 2 /M 2 )4 2 )( A1 ) 

+ c( 4 ) ( g 2 / M 2 )^|M + o(i J ). (i) 



While twist-4 effects are negligible for very large 
Q 2 , at energies of a few GeV they could be rele- 
vant and refine the present QCD predictions. In 
fact, QCD cannot be tested unambiguously un- 
less these contributions are known. 

For twist 4, the Wilson coefficients are known 
at leading order from continuum QCD but 
this does not exclude the possibility to compute 
them also non-perturbatively on the lattice 
The real physical non-perturbative effects are 
contained in the matrix elements of the opera- 
tors in the OPE, and the main problem is the 
lack of quantitative knowledge about these ma- 
trix elements. 



"Talk given by S. Capitani at Lattice 98, Boulder (US). 



2 We do not consider renormalon ambiguities, instanton 
and other possible non-perturbative effects here. 
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For the evaluation of these effects a non triv- 
ial choice of operator bases for higher twist op- 
erators needs to be made The set of all 
possible twist-4 operators is in fact an overcom- 
plete set, and there are different ways to elim- 
inate the redundant operators (using the equa- 
tions of motion) . A "canonical" basis proposed by 
Jaffe and Soldate [|| involves totally symmetric 
and traceless operators with no contracted deriva- 
tives, but even this basis contains more operators 
than needed, as only part of them actually ap- 
pear in the expansion of the T-product at tree 
level. The relevant basis of non-singlet operators 
for minimal spin is^J: 

ijF^F^jpij ipe^PF^F^^ . 

One can also choose a (non-minimal) basis 
in which operators do not contain gluon field 
strengths explicitly, but then contracted deriva- 
tives appear [||. Most likely, the practical choice 
of a basis is to be done case by case, looking at 
the complexities of the mixing patterns, at the 
goodness of the signals and at other practical 
issues. The situation can be further complicated 
by additional mixings due to the lattice. 

Sometimes one can simplify the mixing struc- 
tures. Purely gluonic operators can be at once 
eliminated by including only non-singlet opera- 
tors from the start (as in the basis above), but 
there also exist particular non-trivial combina- 
tions of operators that have protection from sin- 
glet mixing. As an example, the difference be- 
tween the second moments of longitudinal struc- 
ture functions M£(2, Q 2 ) - ^M|(2, Q 2 ) is a pure 
non-singlet quantity also at twist-4 level [Q. 

Of special importance in the twist-4 case are 
the 4-fermion operators. At 1-loop level they do 
not mix with operators of the kind ipFip, and the 
mixing matrix is in this respect triangular. Fur- 
thermore, only 4-fermion operators can transform 
as the flavor 27-plet, and this can be exploited 
to reduce strongly the room for renormalization 

^Flavor structures are not shown. Purely gluonic oper- 
ators like F/j,aDpD^F°^ belong to the singlet sector and 
appear at higher order in the coupling constant. 



mixing, since the flavor 27-plet is not present at 
twist-2 level (where only the singlet and the octet 
can contribute) . Flavor symmetries can be in gen- 
eral very useful: for example the isospin-2 combi- 
nation of the pion structure functions j(| , 

Fi= 2 =F„+ +F„- -2F w o, (2) 

gets contributions only from the 27-plet, and 
therefore cannot mix with twist-2 operators. 

When mixing with operators of leading twist is 
forbidden, one can also avoid renormalon ambi- 
guities completely. But even in cases where this 
mixing is not prohibited, a computation of the 
coefficient functions in a non-perturbative way || 
would get rid of renormalon problems. 

3. PERTURBATIVE 

RENORMALIZATION 

Renormalization factors are needed to relate 
the numbers obtained from lattice simulations 
to the corresponding physical matrix elements: 
O cont {n) = Z{na,g{ii)) ■ O latt {a). In this way 
one can give "continuum" numbers, in the sense 
that a primary result obtained from the lattice 
can be "converted" (through these renormaliza- 
tion factors) to its continuum equivalent. 

Perturbative lattice renormalization is impor- 
tant by itself and as a hint and reference for non- 
perturbative renormalization studies, especially 
when one has to understand lattice mixings, gen- 
erally much more intricate than in the continuum 
case. Mixing patterns on the lattice are in fact 
more transparent when looked at in perturbation 
theory. 

Here we would like to study the perturbative 
renormalization for an important class of twist 4 
operators: the 4-quark operators. We consider 
the symmetrized operators 

°{^} = ^2 ^7{^"V ' i>lv}t A i> ~ traces (3) 

A 

= H ^l{^t A ^ ■ iplv}lst A ip - traces. 

A 

We impose as renormalization conditions that the 
1-loop amputated matrix elements renormalized 
at a reference scale ji are equal to the correspond- 
ing bare tree-level quantities: 
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(q(p)q(p')\G(»)\q(p)q(p')) 



(4) 



Similar results are found for the operator O 



1— loop 

p 2 —p' 2 —^ 2 



Zo (q(p)q(p')\0(a)\q(p)q(p')) 



(q(p)q(p')\0(a)\q(p)q(p')) 



We find it convenient to use forward matrix ele- 
ments (p = p'), the renormalization factors being 
not affected by this choice, and the calculations 
being much simplified. 

The computations have been performed using 
Form codes for the analytic part, and Fortran 
codes for the numerical integration. They have 
also been checked by hand. In principle Fierz 
transformations could be used (both for color and 
Dirac indices), together with charge conjugation 
transformations. However, they are not needed 
when suitable sets of diagrams are evaluated to- 
gether. We find that in this way the computations 
are simpler, and since we can integrate numeri- 
cally each diagram with a good precision in a very 
short time we are not bound to use the 2-fcrmion 
results to which the 4-fcrmion diagrams can be 
reduced by Fierz rearrangements. Furthermore, 
we avoid any problem related to <i-dimensional 
Fierz transformations. 

From our lattice results we can derive the Z fac- 
tors to any continuum scheme. As an example, we 
give here the matching factors between the lattice 
and the MS scheme at a scale fi = 1/a. With ob- 
vious notations, the result for the operator O^^y 
(with fi^i/) is: 

z latt^MS = QA ®tAy X _ g 2 . (5) 

■\{t A ® ^)(-0.001442 + 0.040022 • (1 - a)) 
+ (1 <g> l)(-0.023980 + 0.021345 • (1 - a))}, 

where a is the gauge parameter. To complete the 
mixing matrix we also need the contribution 



rT-latt—tMS 
Z ll 



(6) 



(1 ® 1) • 1 - gl ■ 
•{(1 ® 1) (-0.348170) 

+ (t A <8> ^)(-0.107911 + 0.096053 • (1 - a))\ 



7 latt~^MS 
J tt 



(5) . 

(7) 



■{{t A O t A )(— 0.011619 + 0.040022- (1 - a)) 
+ (1 ® l)(-0.023980 + 0.021345 ■ (1 - a))\ 

z latt^US = (1 g, 1) . 1 _ p 2 . (8) 

•{(1 <g l)(-0.266750) 

+{t A ® t A ){-0. 107911 + 0.096053 • (1 - a))}. 

For these calculations we have used a regulariza- 
tion scheme which involves a totally anticommut- 
ing 75 ({75, 7^} = 0), that is simple to implement 
in computer codes. We are considering perform- 
ing the calculations also in the 't Hooft-Veltman 
scheme. 

The calculation of the renormalization factors 
can also be done non-perturbatively 0], and it 
will be interesting to compare the perturbative 
with the non-perturbative results. Other twist-4 
operators will also be considered, and a choice of 
basis will have to be made, maybe according to 
the physical situations, so that a reliable estimate 
of twist-4 effects for some particular processes can 
be obtained. 
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